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Abstract 

A theoretical  investigation  of  surface  magneto- 
plasmons  in  n-type  InSb  has  been  carried  out  with  the 
effect  of  damping  taken  into  account.  Part  of  the 
investigation  included  surface  magnetoplasmons  inter- 
acting with  surface  optical  phonons.  The  calculated 
dispersion  curves  show  backbending,  in  contrast  to  the 
situation  where  damping  is  neglected.  In  addition, 
there  are  gaps  in  the  dispersion  curves  foi\certain 

frequency  ranges. 

. \ 


* Work  supported  in  part 
under  Grant  No.  N00014 


nV 


: □ 


a a 


by  the  Office  of  Naval  Research 
76-C-0121 


Technical  Report  #78-9 


DSR3im'?VA;  3TE5 


2. 


! 


> 


INTRODUCTION 

Brion  et  have  made  a theoretical  investigation  of 

surface  magnetoplasmons  in  n-type  InSb  with  the  magnetic  field 

parallel  to  the  surface  and  the  wave  vector  perpendicular  to  the 

magnetic  field  direction.  They  obtained  the  interesting  result 

that  a gap  may  appear  in  the  surface  magnetoplasmon  dispersion 

curve;  however,  they  neglected  the  effect  of  damping.  Brion 
( 2^ 

et  al>  have  also  made  a theoretical  investigation  of  the  effects 
of  a magnetic  field  on  interacting  surface  plasmons  and  surface 
optical  phonons,  again  neglecting  damping. 

Experimentally,  damping  is  always  present  due  to  the  scatter- 
ing of  the  free  carriers  by  defects  or  phonons.  In  the  present 
investigation,  we  have  therefore  obtained  the  surface  magnetoplasmon 
dispersion  relation  in  the  presence  of  free  carrier  damping.  The 
interaction  with  surface  optical  phonons  is  also  considered.  Dis- 
persion curves  are  presented  for  n-type  InSb,  using  the  same 

(1  2) 

values  of  the  pertinent  parameters  as  were  used  by  Brion  et  al.  ’ 

The  results  are  compared  with  those  of  Brion  et  al. 


THEORY 


We  will  first  obtain  the  dielectric  tensor  including  damping. 

In  what  follows,  the  static  magnetic  field  5' is  taken  to  be  parallel 
to  the  surface  (which  lies  in  the  xy-plane)  and  pointing  in  the 
y-direction.  In  addition,  we  take  the  magnetoplasmon  wave  vector  to 
be  perpendicular  to  The  geometrical  arrangement  is  shown  in  Fig.  1. 

We  begin  with  the  equation  of  motion 


-e2  - |(v  X 30) 


* dv  * - 
m ^ + m vv 


(1) 
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where  m*  is  the  electron  effective  mass,  v is  the  velocity,  v is 
the  damping  constant,  Eis  the  electric  field,  and  -e  is  the 
electronic  charge.  It  is  assumed  that 


v 


(2) 


and 


E 


-iu)  t 
e 


(3) 


Substituting  these  values  into  Eq.  (1),  we  solve  for  v in  terms 
of  li.  The  conductivity  tensor  a can  then  be  obtained  from  the 
defining  equation 

j = - NQev  = a • E ( 4 ) 

where  N is  the  free  carrier  density  and  j is  the  current 
o 

density. 

The  dielectric  tensor  follows  from  the  equation 

£ - €P  1 + ini  a , (5) 

where  contains  contributions  to  the  dielectric  tensor  other 
than  from  free  carriers  and  1 is  the  3x3  unit  tensor.  In 
component  form,  we  have 


4. 


where 


_ € uj  ”* 

€ (uu)  - ep  - • 


uj  + iv 


XX 


UJ 


(uj+iv)  “ uj 


2 2 ' 


(7a  ) 


i€  u)  ‘ 

6 (cu) 


uu. 


xz 


UJ 


(uj+iv)  2 - uj2 
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- ezx<»>  , 


(7b  ) 


€ uj  i 

e (u)n  - ep  - ” £ • 1 - 

«yyVU);  t u»  UJ+iV 


(7c  ) 


In  these  expressions,  oj  is  the  cyclotron  frequency  given  by 

c 

4c  2 * 

eBQ/m  c,  ujp  is  the  plasma  frequency  given  by  4rrNoe  /m  and 
6 is  the  background  dielectric  constant. 

CD 

Two  values  for  the  quantity  €P  are  considered.  The  situation 
where  the  interaction  between  surface  optical  phonons  and  magneto- 
plasmons  occurs  can  be  handled  in  a simple  manner  by  setting 


( 8) 


where  <uL  and  ujt  are  the  long  wavelength  longitudinal  and  trans- 
verse optical  phonon  frequencies,  respectively.  Damping  of  the 
optical  phonons  is  neglected  because  this  is  small  compared  to 
the  electronic  damping,  for  the  case  considered.  . In  addition, 
the  limiting  case 


5. 


( 9 ) 


is  also  considered.  This  corresponds  to  the  situation  where  the 
effect  of  optical  phonons  is  neglected. 

Consider  next  the  electric  and  magnetic  fields  of  the 
radiation.  The  electric  field  components  will  be  expressed  as 

ik  x - iuit 

Ea(x,t)  - Ea(z)e  , ( 10  ) 

with  similar  forms  for  the  displacement  field  D and  magnetic  field 
3.  Substituting  these  forms  into  the  v x E and  v x 3 Maxwell 
equations  and  eliminating  3,  one  obtains  equations  relating  the 
x-,  y-,  and  z-components  of  E and  D.  To  complete  the  set  of 
equations,  one  needs  the  relations  between  D and  E in  the  crystal 
(z  > 0)  and  in  the  vacuum  (z  < 0);  these  are 


D (z)  =■  E 6 ft(u))Ea(z) 
® p aP  p 

and 

Dtf(z)  = Eff(z)  , 


(11  ) 

(12  ) 


respectively. 

For  surface  magnetoplasmons  for  the  geometry  stated,  we  need 

j 

consider  only  p-polarization,  i.e.,  the  electric  vector  lies  in  the 
sagittal  plane.  For  this  situation,  E (z)  - 0,  and  we  obtain  the 

y 

following  coupled  differential  equations: 


i 


- ^ Vz)  + ikx  33  Ez(z)  - 7 €XX(“-)EX(Z)+  4 €xz(l“)Ez(z) 


(13  ) 


t 4-  E (z)+k2  E < 

x dz  xv  x z 


ikv  £ Ex(2)+k^  Ez(z)  ” 4 €zx(tu)Ex(z)  +S!i2  €zz(u,)Ez(z) 


(14  ) 


If  it  is  assumed  that 


E (z)  = A1e"aZ,  Ez(z)  = A2e"aZ,  z a 0; 


(15  ) 


one  obtains  from  Eqs.  (13)  and  (14)  the  following  expression  for  the 


decay  constant  a: 


2 - k2 
a = k 


-$( 


€xx+  r 


■)  ■ 


(16  ) 


The  constants  A and  A0  are  related  by  the  equation 

A2  " Al(  o 2 / 

\kx  ' / 

c 


(17  ) 


Consider  next  the  fields  in  the  vacuum.  Proceeding  as  before 


and  assuming  that 


Ex(z)  - A3e  ° , Ez(z)  - Dz(z)  - A4e  ° , z<0,  (18  ) 


one  obtains  the  results 


(19  ) 


— 21  a3 

“o  3 


(20  ) 


7. 


By  applying  the  boundary  conditions  of  continuity  of  Ex(z) 

and  D (z)  at  z - 0,  we  obtain  the  surface  raagnetoplasmon  dis- 
z 

persion  relation 


ikx6xz  + or€xx 


k2  - 


2xx 


= 0 . 


(21  ) 


RESULTS  AND  DISCUSSION 

Theoretical  dispersion  curves  have  been  obtained  for  surface 
magnetoplasmons  where  the  effects  of  both  damping  and  retardation 
have  been  included.  The  situation  whe  -"ace  optical  phonons 

interact  with  damped  surface  magnetoplasnw  s also  been  in- 

vestigated. 


Calculations  were  made  using  values  of  parameters  appropriate 


for  n-vype  InSb.  For  the  case  of  surface  magnetoplasmons  alone, 

which  will  be  considered  first,  we  take  €P  * € , i.e.,  we  neglect 

00 

the  contribution  of' surface  optical  phonons.  The  value  used  for 

(3} 

the  background  dielectric  constant  6^  was  15.68v  7 and  the  value 


of  the  ratio  |uJc/cUp|  was  taken  to  b©  0.5.  Two  values  of  the  damping 


parameter  were  considered:  v - 0.01  u>p  and  0.1uup.  For  comparison, 
the  results  of  Brion  et  al^’2^  for  v " 0 are  also  presented. 


In  making  the  calculations,  real  values  of  uu/u>D  were  assumed 
and  complex  values  of  the  wave  vector  kx(“k^  + ikg^  and  the  decay 
constant  + icr were  calculated  using  a digital  computer. 


8. 


The  values  of  uu/uu^  were  taken  in  steps  of  0.01  or  smaller  for  the 

range  0 s uj/uj  s 1.75.  The  computer  program  calculated  the  complex 
P (4) 

values  of  kx  using  Muller's  method  , a process  of  successive 
approximations  from  initial  guesses  (which  need  not  be  externally 
supplied)  . 

Figure  2 shows  the  surface  magnetoplasmon  dispersion  curves 

for  v/uu  = 0.01  for  both  uu  >0  and  ui  <0.  Figure  3 shows  the 
p c c 

corresponding  frequency  dependence  of  the  imaginary  part  of  the 
wave  vector  associated  with  surface  magnetoplasmons . Figures  4 
and  5 give  the  same  type  of  results  for  v/uup  = 0.1.  Figures  2 and  4 

also  show  the  dispersion  curves  for  v = 0 (dash-dash  curves) . 

The  interaction  of  surface  magnetoplasmons  and  surface  optical 
phonons  is  considered  next.  This  interaction  occurs  for  electron 
concentrations  of  £ 10  cm  The  plasma  frequency  mp  was  chosen 
so  that  the  transverse  optical  phoncn  frequency  lies  in  the  gap 
of  the  surface  magnetoplasmon  dispersion  curve  for  v * 0 (see  Fig.  2) 
Figures  6 and  7 show  the  dispersion  curves  for  tu  >0  and  <ju  <0, 

G C 

respectively,  for  v/tup  = 0.1.  As  in  the  preceding  figures,  the 
curves  for  v = 0 are  also  shown.  Figure  8 shows  the  frequency 
dependence  of  the  imaginary  part  of  the  wave  vector  associated 
with  surface  magnetoplasmons  interacting  with  surface  optical  phonons 
In  Fig.  2,  for  v ■ 0,  we  see  that  for  the  two  signs  of  ujc  (or, 
equivalently,  the  two  signs  of  k^)  , the  dispersion  curves  differ 
for  most  of  the  frequency  range  considered.  Consider  first  the 
case  for  id  >0.  The  dispersion  curve  consists  of  two  branches 

G 

with  a gap  between  them,  i.e.,  the  single  branch  that  would  occur 
for  zero  magnetic  field  (not  shown)  is  split  into  two  branches  when 


9. 


a sufficiently  large  magnetic  field  is  applied.  The  lower  branch^ ^ 
of  the  surface  magnetoplasmon  dispersion  curve  terminates  when  it 
intersects  the  dispersion  curve  for  bulk  magnetoplasmons  defined 
by  a'  =0  where  a1  is  the  value  of  a when  v =0.  The  upper 
branch^ ^ begins  on  the  line  defined  by  € =*  0,  and  approaches 

aA 

asymptotically  the  frequency  defined  by  the  equation 
1 + £ - i€  “0,  this  being  the  frequency  for  unretarded  surface 

magnetoplasmons  found  by  Pakhomov  and  Stepanov ^ . Brion  et  al^^ 
have  shown  that  for  the  gap  to  occur,  the  dielectric  constant  6^ 

2 2 A 

must  satisfy  the  inequality  6 > (u)  + uu„)2/ou  . 

» c p c 

Consider  next  the  case  for  uj  <0.  Again  for  non-zero  magnetic 

c 

field  there  are  two  branches  in  the  surface  magnetoplasmon  dispersion 

curve  with  a gap  separating  them.  The  lower  branch^^  approaches 

the  asymptotic  frequency  value  given  by  the  equation 

1 + € +i€_“0.  The  upper  branch^  starts  at  the  light  line 

at  the  point  where  € - 1,  and  terminates  when  it  intersects  the 

upper  bulk  magnetoplasmon  dispersion  curve  a'  =0.  Figure  2 shows  that 

the  addition  of  damping  (v,=0.01u)p)  causes  changes,  relative  to  the 

v “ 0 situation,  for  both  signs  of  w . For  uu  >0,  there  are  three 

c c 

branches.  The  lower  branch  coincides  with  the  corresponding  v = 0 
branch  until  the  former  terminates  at  ”0,  just  before  inter- 
secting the  line  a'-  0 (a^  is  nearly  equal  to  since  the  damping 
v is  small) . Then  there  is  a gap  followed  by  the  second  branch 
which  begins  at  a1  ■ 0,  just  to  the  left  of  the  light  line,  and 
continues  upward  until  it  intersects  the  light  line,  at  which  point 
it  terminates.  There  is  then  another  gap  followed  by  the  third 
branch  which  begins  at  the  light  line  near  ^xx  “ 1 and  essentially 


10. 


follows  the  corresponding  v “ 0 curve  as  it  approaches  the  line 


1 + € - i€  - 0 . 

XX  xz 

For  uj  <0,  there  are  two  branches.  The  lower  branch  is 
c 

essentially  coincident  with  the  v = 0 curve  until  it  almost  reaches 

the  asymptotic  line  l+€  = i€  _ “ 0.  Then  it  quickly  bends 

XX  xz 

back,  crosses  the  light  line,  and  terminates  at  a point  where  k0  “ 0 , 
i.e.,  where  the  imaginary  part  of  the  surface  polariton  wave  vector 
vanishes.  A large  gap  occurs,  and  then  the  second  branch  begins 
at  the  light  line  at  k„  = 0,  near  6 =1.  It  initially  follows  the 

v “ 0 curve,  then  backbends  across  the  light  line,  doubles  back 
slightly,  and  subsequently  proceeds  nearly  parallel  to  the  light 
line.  The  backbending  is  analogous  to  the  situation  that  occurs 
for  Ge(6). 

Figure  3 shows  the  frequency  dependence  of  the  imaginary  part 
of  the  wave  vector  for  both  signs  of  ui  . Note  by  comparing  Figs.  2 
and  3 that  when  the  k^  curves  reverse  direction,  the  k2  curves 
also  reverse  direction.  In  a sense,  the  k2  curves  track  the  k^ 
curves.  This  type  of  behavior  has  also  been  noted  before v . 

Figure  4 shows  the  effect  on  the  surface  magnetoplasmons  of 
increasing  the  damping  by  one  order  of  magnitude,  i.e.,  increasing 
the  damping  to  v ■ O.lui  . Note  that,  as  compared  to  the  result  for 

Mr 

v “ O.Oluip  in  Fig.  2,  the  curve  for  v ” 0.1ujp  (<flc>0)  in  Fig.  4 fills 
in  the  gaps  between  the  first  and  second  and  between  the  second  and 
third  branches,  and  then  terminates  to  the  left  of  the  light  line 
at  - 0.  The  lower  branch  of  the  ujc  < 0 curve  moves  to  the  right 
of  the  light  line  to  a much  lesser  extent  for  v/u)p  “0.1  (Fig.  4) 


A 


. * \i‘,  #4 
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than  for  v =*  0.01  uu  (Fig.  2)  before  backbending  and  terminating  at 

P 

k - o.  The  upper  branch  of  the  uu  <0  curve  for  v **  0.1  uu  begins 
2 c p 

at  the  light  line  at  k2  = 0 near  €xx  “ 0,  moves  upward  near  the 

light  line,  and  then  backbends  slightly.  Subsequently,  it  moves 

essentially  parallel  to  the  light  line  in  the  same  manner  as  for 

v = 0.01  uu  (Fig.  2)  . 

P 

The  frequency  dependence  of  the  imaginary  part  of  the  wave 
vector  shown  in  Fig.  5 for  v 
behavior  of  the  frequency  dependence  of  the  real  part  of  the  wave 
vector  shown  in  Fig.  4. 

Next,  consider  the  interaction  of  surface  magnetoplasmons 

( 2) 

and  surface  optical  phonons  in  n-type  InSb.  Brion  et  al.  found 
that  this  interaction,  for  the  case  of  v = 0,  leads  to  interesting 
results  if  the  unretarded  surface  optical  phonon  frequency  uug. 


0.1  uUp  essentially  tracks  the 


given  by 


uu 


( 2 2 V 

2 *“L  ' “T 

u<T  + / 

l i + i/O 


(22  ) 


is  located  in  the  gap  region  for  surface  magnetoplasmons  alone. 

r 

Figure  6 shows  the  dispersion  curves  (uu  > 0)  that  result 

c 

when  the  magnetoplasmon-optical  phonon  interaction  is  included. 
For  v - 0,  the  interaction  has  the  effect  of  doubling  the  number 

of  branches  (compare  Fig.  6 with  Fig.  4).  In  Fig.  6,  the  lower 

* 

branch  starts  at  zero  and  moves  to  the  right  of  the  light  line 
until  it  intersects  a bulk  dispersion  curve  (one  of  the  roots 
of  or'  m 0)  . The  second  branch  starts  at  the  light  line  at  one 


fi  4ft. 
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of  the  roots  of  ^ “0  and  then  flattens  out  asymptotically  as 

XX 

indicated  in  the  figure.  The  third  branch  starts  at  the  light 
line  at  (^(excluding , however,  the  point  uj  ” uu^.)  and  terminates 
when  it  intersects  a bulk  dispersion  curved*7  “ 0).  Finally,  the 
fourth  branch  begins  at  the  light  line  at  € - 0 and  then  flattens 

AA 

out  asymptotically  as  indicated. 

The  inclusion  of  damping  (v  = 0.1uup),  as  is  now  to  be  ex- 
pected, causes  the  dispersion  curves  to  backbend  across  the  light 
line.  In  addition,  there  are  now  only  two  branches  instead  of 
four.  The  lower  branch  at  first  coincides  with  the  v =■  0 branch, 
backbends  in  the  vicinity  of  a'  ” 0,  crosses  the  light  line,  has 
a dimple  which  carries  it  just  to  the  right  of  the  light  line 
at  a frequency  just  below  ujt,  and  then  terminates  at  u)  =*  uuT 
(also  at  this  point,  kg  “0).  The  next  branch  starts  at  the 
light  line  at  w - Wj  and  initially  follows  the  v = 0 curve. 

Near  a ' =0,  however,  this  branch  backbends,  crosses  the  light 
line  to  the  left,  recrosses  the  light  line  to  the  right,  and 
finally  backbends  to  the  left  of  the  light  line  terminating 
at  kg  - 0. 

Figure  7 shows  the  surface  magnetoplasmon-surf ace  optical 
phonon  dispersion  curves  for  u>c  < 0.  As  for  the  case  just  dis- 
cussed Cu)  > 0)  , the  surface  magnetoplasmon-surf  ace  optical 
c 

phonon  interaction  results  in  doubling  the  number  of  branches 
for  v - 0 (compare  Fig.  7 with  Fig.  4).  In  addition,  the 
branches  behave  in  a manner  very  similar  to  that  where  cu  > 0. 

v 

The  addition  of  damping  for  cu  <0  (Fig.  6)  leaves  the 

c 

number  of  branches  at  four.  The  first  branch  initially  tracks 

t 

the  v “ 0 case;  it  then  backbends  and  terminates  at  kg  ■ 0. 

I 
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The  second  branch  is  short,  starting  just  to  the  left  of  the 
light  line  where  k2  - 0,  showing  a dimple  which  carries  it  just 
to  the  right  of  the  light  line  at  a frequency  just  below  uuT, 
and  terminating  at  uj  ■ ujt,  where  also  k2  = 0.  The  third 
branch  starts  at  m » tuT,  moves  slightly  to  the  right  of 

the  light  line,  backbends  and  then  terminates  at  k2  - 0.  Finally, 
the  fourth  branch  begins  at  the  light  line  at  k2  - 0,  and  moves 
essentially  along  the  light  line,  but  slightly  to  the  right.  It 
subsequently  backbends  to  the  left  of,  and  then  moves  parallel  to, 
the  light  line. 

Figure  8 shows  the  corresponding  frequency  dependence  of  the 
imaginary  part  of  the  wave  vector.  The  behavior  is  similar  to 
that  noted  previously. 

Although  we  have  not  carried  out  an  exhaustive  investigation 

of  the  effects  of  damping  on  surface  magnetoplasmon  dispersion 

curves,  our  results  demonstrate  that  damping  produces  backbending 

C 6 7) 

of  the  type  found  in  previous  investigations v * ' when  the  fre- 
quency is  taken  to  be  real  and  the  wave  vector  is  taken  to  be 
complex.  As  a consequence  of  the  backbending,  portions  of  the 

dispersion  curves  lie  to  the  left  of  the  light  line  ui  ■ ck,  and  the 

(8) 

corresponding  modes  are  known  as  Zenneck  modes. 

Another  effect  produced  by  damping  is  a tendency  to  reduce 

or  eliminate  gaps  in  the  dispersion  curves  as  the  damping  is 

increased.  This  is  clearly  illustrated  in  Fig.  4 for  v “ O.lujp 

where  the  gap  in  the  dispersion  curve  for  uj  >0  has  been 

c 

completely  eliminated. 
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Some  comment  on  the  relationship  of  our  results  to  the 

(Q) 

experimental  work  of  Hartstein  and  Burstein  is  in  order. 

In  their  determination  of  the  surface  magnetoplasmon  dispersion 

curves  for  n-type  InSb  using  attenuated  total  reflection,  they 

observed  no  backbending.  This  result  may  be  understood  in 

(7) 

terms  of  the  discussion  of  Kovener  et  al:  who  point  out  that 

dispersion  curves  determined  by  varying  the  incident  angle 
at  fixed  uj  should  show  backbending,  while  those  determined  by 
varying  yj  at  fixed  incident  angle  should  not  show  backbending. 

Since  Hartstein  and  Burstein  carried  out  their  experiments  by 
varying  uu  at  fixed  incident  angle,  their  procedure  falls  into 
the  second  category  just  mentioned,  and  it  is  reasonable  that 
they  did  not  observe  backbending.  It  would  be  interesting,  however, 
to  carry  out  experiments  to  reveal  the  backbending  by  varying  the 
incident  angle  at  fixed  uu . 
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FIGURE  CAPTIONS 

1.  The  coordinate  system  and  specimen  geometry. 

2.  Dispersion  curves  for  surface  magnetoplasmons  in  n-type  InSb 

<v/»  " 0.01).  The  dash  - dash  curves  are  the  dispersion 

curves  for  v/u)  - 0.0. 

P 

3.  The  frequency  dependence  of  the  imaginary  part  of  the  wave 
vector  associated  with  surface  magnetoplasmons  in  n-type  InSb 
(v/oup  - 0.01). 

4.  Dispersion  curves  for  surface  magnetoplasmons  in  n-type  InSb 

(v/ujp  - 0.1).  The  dash  - dash  curves  are  the  dispersion 

curves  for  v/ui  - 0.0. 

P 

5.  The  frequency  dependence  of  the  imaginary  part  of  the  wave 

vector  associated  with  surface  magnetoplasmons  in  n-type 

InSb  (v/cu  " 0.1). 

P 

6.  Dispersion  curves  for  surface  magnetoplasmons  In  n-type  InSb 
where  the  interaction  with  surface  optical  phonons  is  included 
(uu  >0).  The  dash  - dash  curves  are  the  dispersion  curves  for 

v/tUp  ■ 0.0. 

7.  Dispersion  curves  for  surface  magnetoplasmons  in  n-type  InSb, 
where  the  Interaction  with  surface  optical  phonons  is  is  included 
(<«  <0).  The  dash  - dash  curves  are  the  dispersion  curves  for 

C 

v/«)p  "0.0. 

8.  The  frequency  dependence  of  the  imaginary  part  of  the  wave 
vector  associated  with  surface  magnetoplasmons  in  n-type  InSb 
where  the  interaction  with  surface  optical  phonons  is  included. 
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